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We compute the entanglement entropy for some quantum field theories on de Sitter 
space. We consider a superhorizon size spherical surface that divides the spatial slice into 
two regions, with the field theory in the standard vacuum state. First, we study a free 
massive scalar field. Then, we consider a strongly coupled field theory with a gravity dual, 
computing the entanglement using the gravity solution. In even dimensions, the interesting 
piece of the entanglement entropy is proportional to the number of e-foldings that elapsed 
since the spherical region was inside the horizon. In odd dimensions it is contained in a 
certain finite piece. In both cases the entanglement captures the long range correlations 
produced by the expansion. 



1. Introduction 

Entanglement entropy is a useful tool to characterize states with long range quantum 
order in condensed matter physics (see and references therein). It is also useful in 
quantum field theory to characterize the nature of the long range correlations that we have 
in the vacuum (see e.g. 0,^0 and references therein). 

We study the entanglement entropy for quantum field theories in de Sitter space. 
We choose the standard vacuum state (the Euclidean, Hartle-Hawking/Bunch- 

Davies/Chernikov-Tagirov vacuum). We do not include dynamical gravity. In particular, 
the entropy we compute should not be confused with the gravitational de Sitter entropy. 

Our motivation is to quantify the degree of superhorizon correlations that are gener- 
ated by the cosmological expansion. 

We consider a spherical surface that divides the spatial slice into the interior and 
exterior. We compute the entanglement entropy by tracing over the exterior. We take 
the size of this sphere, i?, to be much bigger than the de Sitter radius, R ^ Rds = , 
where H is Bubble's constant. Of course, for R <^ Rds we expect the same result as in 
fiat space, li R = Rds-, then we would have the usual thermal density matrix in the static 
patch and its associated entropy^. As usual, the entanglement entropy has a UV divergent 
contribution which we ignore, since it comes from local physics. For very large spheres, 
and in four dimensions, the finite piece has a term that goes like the area of the sphere and 
one that goes like the logarithm of the area. We focus on the coefficient of the logarithmic 
piece. In odd spacetime dimensions there are finite terms that go like positive powers of 
the area and a constant term. We then focus on the constant term. 

We first calculate the entanglement entropy for a free massive scalar field. To de- 
termine it, one needs to find the density matrix from tracing out the degrees of freedom 
outside of the surface. When the spherical surface is taken all the way to the boundary of 
de Sitter space the problem develops an S'0(1, 3) symmetry. This symmetry is very helpful 
for computing the density matrix and the associated entropy. Since we have the density 
matrix, it is also easy to compute the Renyi entropies. 

We then study the entanglement entropy of field theories with a gravity dual. When 



the dual is known, we use the proposal of [10,11] to calculate the entropy. It boils down to 



an extremal area problem. The answer for the entanglement entropy depends drastically 



^ This can be regarded as a (UV divergent) 0{G%) correction to the gravitational entropy of 
de Sitter space S. 
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on the properties of the gravity dual. In particular, if the gravity dual has a hyperbolic 
Friedman-Robertson- Walker spacetime inside, then there is a non-zero contribution at 
order A^^ for the "interesting" piece of the entanglement entropy. Otherwise, the order A^^ 
contribution vanishes. 

This provides some further hints that the FRW region is indeed somehow contained 
in the field theory in de Sitter space [|T^. More precisely, it should be contained in the 
superhorizon correlations of colored fieldsil. 

The paper is organized as follows. In section 2, we discuss general features of entan- 
glement entropy in de Sitter. In section 3, we consider a free scalar field and compute 
its entanglement entropy. In section 4, we write holographic duals of field theories in de 
Sitter, and compute the entropy of spherical surfaces in these theories. We end with a 
discussion. Some more technical details are presented in the appendices. 



2. General features of entanglement entropy in de Sitter 

Entanglement entropy is defined as follows 0. At some given time slice, we consider 
a closed surface E which separates the slice into a region inside the surface and a region 
outside. In a local quantum field theory we expect to have an approximate decomposition 
of the Hilbert space into H = Hin x Hout where Hin contains modes localized inside 
the surface and Hout modes localized outside. One can then define a density matrix 
Pin = Tr/fo„t IV') ('01 obtained by tracing over the outside Hilbert space. The entanglement 
entropy is the von Neumann entropy obtained from this density matrix: 

S = -Tr Pin log Pin (2.1) 

2.1. Four dimensions 

We consider de Sitter space in the fiat slicing 

ds^ = , i-drf + dx\ + dxl + dxl) (2.2) 



A holographic calculation of the entanglement entropy associated to a quantum quench is 



presented in |13]. A quantum quench is the sudden perturbation of a pure state. The subsequent 
relaxation back to equilibrium can be understood in terms of the entanglement entropy of the 
quenched region. There, one has a contribution to the (time dependent) entropy coming from the 
region behind the horizon of the holographic dual. 
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where H is the Hubble scale and ry is conformal time. We consider surfaces that sit at 
constant ry slices. We consider a free, minimally coupled, scalar field of mass m in the 
usual vacuum state 

As in any quantum field theory, the entanglement entropy is UV divergent 

S = •S'uv-divergent + 5'uV-finite (2.3) 

The UV divergencies are due to local effects and have the the form 

5uV-divergent = Cl ^ + log(ei?) (C2 + C3 Am^ + C4AH^) (2.4) 

where e is the UV cutoff. The first term is the well known area contribution to the 
entropy coming from entanglement of particles close to the surface considered. The 
logarithmic terms involving C2 and C3 also arise in flat space. Finally, the last term involves 
the curvature of the bulk spacei. All these UV divergent terms arise from local effects and 
their coefficients are the same as what we would have obtained in flat space. We have 
included if as a scale inside the logarithm. This is just an arbitrary definition, we could 



also have used m when m is non-zero. 

Our focus is on the UV finite terms that contain information about the long range 
correlations of the quantum state in de Sitter space. The entropy is invariant under the 
isometries of dS. This is true for both pieces in ( f2.3| ). In addition, we expect that the 
long distance part of the state becomes time independent. More precisely, the long range 
entanglement was established when these distances were subhorizon size. Once they moved 
outside the horizon we do not expect to be able to modify this entanglement by subsequent 
evolution. Thus, we expect that the long range part of the entanglement entropy should 
be constant as we go to late times. So, if we fix a surface in comoving x coordinates in 
( p.2| ), and we keep this surface fixed as we move to late times, 77 — )■ 0, then we naively 



^ In de Sitter there is only one curvature scale, but in general we could write terms as 

S\ogeH = J [aR^^pan'^n^njnj + bR^j.uu'^n'i +cR + dKf" Ki^^ + eK^^Ki^ ^ ) (2.5) 

where K are the extrinsic curvatures and i,j label the two normal directions and z/, • • • are 
spacetime indices The extrinsic curvatures also contribute to C2 in ( [2.41 ). One could also write a 
term that depends on the intrinsic curvature of the surface, R-e, but the Gauss-Codazzi relations 



can be used to relate it to the other terms in (2.5). 
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expect that the entanglement should be constant. This expectation is not quite right 
because new modes are coming in at late times. However, all these modes only give rise to 
entanglement at short distances in comoving coordinates. The effects of this entanglement 
could be written in a local fashion. 

In conclusion, we expect that the UV-finite piece of the entropy is given by 

-Suv-Finite = csAif^ + — \og{AH^) + finite = + ce log?? + finite (2.6) 

2 rj^ 

where A is the proper area of the surface and A^ is the area in comoving coordinates 
(^ = jj^^a )■ The finite piece is a bit ambiguous due to the presence of the logarithmic 
term. 

The coefficient of the logarithmic term, cg, contains information about the long range 
entanglement of the state. This term looks similar to the UV divergent logarithmic term 
in ( |2.4| ), but they should not be confused with each other. If we had a conformal field 
theory in de Sitter they would be equal. However, in a non-conformal theory they are 
not equal (cg 7^ C2). For general surfaces, the coefficient of the logarithm will depend on 
two combinations of the extrinsic curvature of the surface in comoving coordinates. For 
simplicity we consider a sphere here^. This general form of the entropy, ( |2.6| ), will be 
confirmed by our explicit computations below. 

We define the "interesting" part of the entropy to be the coefficient of the logarithm, 
S'intr = ce. The t/F-finite area term, with coefficient C5, though physically interesting, is 
not easily calculable with our method. It receives contributions from the entanglement at 
distances of a few Hubble radii from the entangling surface. It would be nice to find a 
way to isolate this contribution and compute C5 exactly. We could only do that in the case 
where the theory has a gravity dual. 

2.2. Three dimensions 

For three dimensional de Sitter space we can have a similar discussion. 

S = di h 5'uv-finite 

A, ^2.7) 

'S'uV-finite = d2AH + ^3 = ^2 1" <^3 

?7 



^ It is enough to do the computation for another surface, say a cyhnder, to determine the 
second coefficient and have a result that is vahd for general surfaces |15]. In other words, for a 



general surface we have ce = /i / KabKab + f2 J {Kaa)^ where /i, /2 are some constants and Kab 
is the extrinsic curvature of the surface within the spatial slice. 
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Here there is no logarithmic term. The interesting term is which is the finite piece. So 
we define S'intr = d^. 

A similar discussion exists in all other dimensions. For even spacetime dimensions 
the interesting term is the logarithmic one and for odd dimensions it is the constant. One 
can isolate these interesting terms by taking appropriate derivatives with respect to the 
physical done in in a similar contexts. 

Note that we are considering quantum fields in a fixed spacetime. We have no gravity. 
And we are making no contact with the gravitational de Sitter entropy which is the area 
of the horizon in Planck units. 



3. Entanglement entropy for a free massive scalar field in de Sitter 

Here we compute the entropy of a free massive scalar field for a spherical entangling 
surface. 

3.1. Setup of the problem 

Consider, in flat coordinates, a spherical surface S'^ defined by +0:3 = R^. We 

consider Rc ^ r]. This means that the surface is much bigger than the horizon. 

If we could neglect the ry dependent terms, we can take the limit 77 — > 0, keeping Rc 
fixed. This then becomes a surface on the boundary. This surface is left invariant by an 
SO{l, 3) subgroup of the 5*0(1, 4) de Sitter isometry group. We expect that the coefficient 
of the logarithmic term that we discussed above is also invariant under this group. It 
is therefore convenient to choose a coordinate system where 5'0(1,3) is realized more 
manifestly. This is done in two steps. First we can consider de Sitter in global coordinates, 
where the equal time slices are three- spheres. Then we can choose the entangling surface 
to be the two-sphere equator of the three-sphere. In fact, at 77 = 0, we can certainly map 
any two sphere on the boundary of de Sitter to the equator of by a de Sitter isometry. 
Finally, to regularize this problem we can then move back the two sphere to a very late 
fixed global time surface. 

We can then choose a coordinate system where the S0(l,3) symmetry is realized 
geometrically in a simple way. Namely, this S0(l,3) is the symmetry group acting on 
hyperbolic slices in some coordinate system that we describe below. 



See formula (1.1) of |16]. 
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(a) (b) (c) 

Fig. 1: Setup of the problem: (a) We consider a spliere with radius much greater 
than the horizon size, at late conformal time rj, in flat slices, (b) This problem can 
be mapped to half of a 3-sphere , also with boundary S"^ , but now the equator, at 
late global time tb ■ (c) We can also describe this problem using hyperbolic slices. 
The interior of the sphere maps to the "left" (L) hyperbolic slice. The Penrose 
diagrams for all situations are depicted below the geometric sketches. 



3.2. Wavefunctions of free fields in hyperbolic slices and the Euclidean vacuum 



The hyperbolic/open slicing of de Sitter space was studied in detail in [^|T8[. It 



can be obtained by analytic continuation of the sphere 5"* metric, sliced by 5'"^s. The 
S"^ is described in embedding coordinates by Xf + ... + X| = H~^. The coordinates are 
parametrized by angles in the following way: 

= coste cos pe, X4^ = H~^smTE: Xi,2,3 = cosr^; sinp^;ni,2,3 (3.1) 

where rii are the components of a unit vector in R^. The metric in Euclidean signature is 
given by: 

ds% = H'^idrl + cos^ TEidp^ + sin^ Pe dQl)) (3.2) 
We analytically continue X5 — )■ iXo. Then the Lorentzian manifold is divided in three 
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parts, related to the Euclidean coordinates by: 

^ . I = f - iti? tR>0 
PE = -irR rR>0 



C: r^"^? ■ -"/2<^^<-/2 (3.3) 

TE = -^+itL tL>0 

PE = -irL rL>0 



The metric in each region is given by: 

dsR = H~'^{-dt]^ + sinh^ tR{dr\ + sinh^ ^^^^02)) 

dsl; = H-^{dt'^c + cos^ tc{-drl + cosh^ rcdfll)) (3.4) 

dsl = H-^{-dtl + sinh^ tLidrl + sinh^ VLdnl)) 

We now consider a minimally coupled^ massive scalar field in ^154, with action given by 
S = ^ f \/—g{—(V4>)'^ — 'ni?(fP'). The equations of motion for the mode functions in the R 
or L regions are 



smh"^ t- n-L?f3 + - - ly 



u{t,r,n) = (3.5) 



.sinh^t^t dt sinh^t " 

Where L^s is the Laplacian in the unit hyperboloid, and the parameter u is 
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4~lP 



(3.6) 



When = \ (or ^ = 2) we have a conformally coupled massless scalar. In this case 
we should recover the fiat space answer for the entanglement entropy, since de Sitter is 

2 

conformally fiat. We will consider first situations where ^ > 2, so that 0<z^<l/2or 
V imaginary. The minimally coupled massless case corresponds to v = 3/2. We will later 

2 

comment on the low mass region, ^ <2orl/2<z^<3/2. 

The wavefunctions are labeled by quantum numbers corresponding to the Casimir on 
and angular momentum on 5"^: 

Uplm ~ ^■^^^ Xp(t)yp;m(r, O2) , -Ij^sYplm = (1 + p'^)Yplm (3.7) 



^ If we had a coupling to the scalar curvature £,R4>'^ , we can simply shift the mass rnlff 
m? + 6^H^ and consider the minimally coupled one. 
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The Ypim are eigenfunctions on the hyperboloid, analogous to the standard spherical har- 
monics. Their expressions can be found in fT^j . 

The time dependence (other than the 1 / sinh t factor) is contained in the functions 
Xp{i)- The equation of motion (|3.5|) is a Legendre equation and the solutions are given in 
terms of Legendre functions Pa{x). In order to pick the "positive frequency" wavefunctions 
corresponding to the Euclidean vacuum we need to demand that they are analytic when 
they are continued to the lower hemisphere. These wavefunctions have support on both 
the Left and Right regions. This gives [111 



Xp, 

' 2^i^ l r(. + zp + i/2) ''"-V2(^"^h^^) - fo73-^^."-i/2(^o^ht,) 

(3.8) 

The index a can take the values ±1. For each a the top line gives the function on the 
R hyperboloid and the bottom line gives the value of the function on the L hyperboloid. 
There are two solutions (two values of a) because we started from two hyperboloids. 
The field operator is written in terms of these mode functions as 

^(.,^/* 5: (3.9) 

To trace out the degrees of freedom in, say, the R space, we change basis to functions 
that have support on either the R or L regions. It does not matter which functions we 
choose to describe the Hilbert space. The crucial simplification of this coordinate system 
is that the entangling surface, when taken to the de Sitter boundary, preserves all the 
isometries of the slices. This implies that the entanglement is diagonal in the p,l,m 
indices since these are all eigenvalues of some symmetry generator. Thus, to compute this 
entanglement we only need to look at the analytic properties of ( |3.8| ) for each value of p. 

Let us first consider the case that u is real. For the R region we take basis functions 
equal to the Legendre functions P^^_-^^2i'^osh.tii) and P~^^^^{cosh.tii), and zero in the L 
region. These are the positive and negative frequency wavefunctions in the R region. 
We do the same in the L region. These should be properly normalized with respect to 
the Klein-Gordon norm, which would yield a normalization factor Np. We can write the 
original mode functions, ( |3.8| ), in terms of these new ones in matricial form: 



m\p^n-^ 



The capital indices (/, J) run from 1 to 4, as we are grouping both the cind Xa- The 
coefficients a and (3 are simply the terms multiplying the corresponding P functions in 
( P?8|) , see appendix A for their explicit values. As the field operator should be the same 
under this change of basis, then it follows that: 

<t> = o.,x' = bjP^N-^ ^aj = b,(M-^)'j 

q — R^L 

Here = (a^r, aj.), = {pL^R, h^^ i?)' ^'^'^ ~ {Pl,r, Pl,r)- M is a 2 x 2 matrix whose 
elements are 2x2 matrices. The expression for is the definition of 5, 7, etc. The 

vacuum is defined so that a^Y^) = 0. We want to write in terms of the hn^L oscillators 
and the vacua associated to each of these oscillators, 6i?|-R) = and 6l|L) = 0. As we are 
dealing with free fields, their Gaussian structure suggests the ansatz 

1^) = e^^^^^=n,L'^^^^>]\R)\L) (3.12) 
and one can solve for rriij demanding that acr\^) = 0. This gives 

^ijlja + = = -5io-(7~"^)crj (3.13) 

Using the expressions in ( |3.8| ) (see appendix A) we find for m: 

^ y2e-P" / cosnu isinhp7r\ 

v/cosh27rp + cos27rz^ \^zsinhp7r costxu J 

Where 6 is an unimportant phase factor, which can be absorbed in the definition of the 
oscillators. In rriij the normalization factors Np drop out, so they never need to be 
computed. 

The expression ( |3.12| ), with ( |3.14D , needs to be simplified more before we can easily 
trace out the R degrees of freedom. We would like to introduce new oscillators cl and cr 
(and their adjoints) so that the original state \1/ has the form 

l^r) =e'^^fi^^|i?)'|L)' (3.15) 

where \R)'\L)' are annihilated by Ci?, cl- The details on the transformation are in ap- 
pendix A. Here we state the result. The 6's and c's are related by: 

cr = ubn + vb^j^ 

t , ,2 , ,2 (3.16) 
CL=ubL + vb[^, \u\ — \v\ =1 
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Requiring that Cfl|\l/) = 7c|^|\l/) and cl|\1/) = 7c]^|^') imposes constraints on u and v. The 
system of equations has a solution with 7 given by 

7 = i = (3.17) 

V cosh 2np + cos 2nh' + y/ cosh 27rp + cos 2nh' + 2 

We have considered the case of < z/ < 1/2. For z/ imaginary, ( |3.17| ) is analytic under 



the substitution v — )• iv, which corresponds to substituting cos 2niy — )• cosh 2nih', so ( |3.17| ) 
is also valid for this range of masses. One can check directly, by redoing all the steps in 
the above derivation, that the same final answer is obtained if we had assumed that u was 
purely imaginary. 

3.3. The density matrix 

The full vacuum state is the product of the vacuum state for each oscillator. Each 
oscillator is labelled by p, /, m. For each oscillator we can write the vacuum state as in 
( p.l5|) . Expanding (|3.15|) and tracing over the right Hilbert space we get 

00 

Pp,i,m = Tr//^(|^')(*|) oc ^ |7pp|n;p,/,m)(n;p,/,m| (3.18) 

n=0 

So, for given quantum numbers, the density matrix is diagonal. It takes the form pl{p) = 
(1 — |7pP)diag(l, |7pP, |7p|'*, ■ ■ ■)' normalized to Trpi, = 1. The full density matrix is simply 
the product of the density matrix for each value of p, /, m. This reflects the fact that there 
is no entanglement among states with different 5'0(1,3) quantum numbers. The density 



matrix for the conformally coupled case was computed before in |T^ . 

Here, one can write the resulting density matrix as = 6"^^"="* with Hent called the 
entanglement hamiltonian. Here it seems natural to choose (3 = 2n as the inverse temper- 
ature of dS. Because the density matrix is diagonal, the entanglement Hamiltonian should 
be that of a gas of free particles, with the energy of each excitation a function of the 
Casimir and the mass of the scalar field. This does not appear to be related to any ordinary 
dynamical Hamiltonian in de Sitter. In other words, take pL oc diag(l, |7pP, |7p|^, ■■■) then 
the entanglement Hamiltonian for each particles is Hp = Epcj^Cp, with Ep = log |7pp. 
For the conformally coupled scalar then Ep = p and we have the entropy of a free gas in 
H^. In other words, in the conformal case the entanglement Hamiltonian coincides with 



the Hamiltonian of the field theory on i? x H [20.21 
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3.4- Computing the Entropy 



With the density matrix ( p.l8|) we can calculate the entropy associated to each par- 



ticular set of 5(9(1,3) quantum numbers 



Sip, u) = -TipM \ogpM = - log(l - |7p|') - , log |7p|' (3.19) 



1-h 



The final entropy is then computed by summing ( |3.19| ) over all the states. This sum 
translates into an integral over p and a volume integral over the hyperboloid. In other 
words, we use the density of states on the hyperboloid: 

S{^) = Vh^ j dpV{p)S{p, jy) (3.20) 

The density of states for radial functions on the hyperboloid is known for any dimensions 

„ 2 

|22} j. For example, for H , T>{p) = Here Vh^ is the volume of the hyperboloid. This 
is of course infinite. This infinity is arising because we are taking the entangling surface 
all the way to = 0. We can regularize the volume with a large radial cutoff in . This 
should roughly correspond to putting the entangling surface at a finite time. Since we are 
only interested in the coefficient of the logarithm, the precise way we do the cutoff at large 
volumes should not matter. The volume of a unit size for radius less that is given 
by 



= j rfr sinh^ r ~ 47r - (3.21) 

The first term goes like the area of the entangling surface. The second one involves the 
logarithm of this area. We can also identify — t- — log 77. This can be understood more 
precisely as follows. If we fix a large t^ and we go to large r^, then we see from (|3.1|)(|3.3| ) 
that the corresponding surface would be at an 77 oc e~'~^, for large r^. Thus, we can 
confidently extract the coefficient cq in ( p.6|) . For such purposes we can define Vf/s^eg = 27r. 
The leading area term , proportional to e^^" depends on the details of the matching of this 
IR cutoff to the proper UV cutoff. These details can change its coefficient. 
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Fig. 2: Plot of the entropy Sintr / Sintr,v=i/2 of the free scalar field, normalized to 
the conformally coupled scalar, versus its mass parameter squared. The minimally 
coupled massless case corresponds to = 9/4, the conformally coupled scalar to 
= 1/4 and for large mass (negative v"^) the entropy has a decaying exponential 
behavior. 



Thus, the final answer for the logarithmic term of the entanglement entropy is 



S = cq log T] + other terms 

"OO 



intr 



C6 = - / dpp S{p, u) 

TT ./o 



(3.22) 



with S{p, u) given in (|3.19|) , (|3.17|) . This is plotted in fig. 2 



3.5. Extension to general dimensions 

These results can be easily extended to a real massive scalar field in any number 
of dimensions D. Again we have hyperbolic H^~^ slices and the decomposition of the 
time dependent part of the wavefunctions is identical, provided that we replace ly by the 
corresponding expression in D dimensions 



(D-l) 



m 



(3.23) 
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Then the whole computation is identical and we get exactly the same function S{p, u) 
for each mode. The final result involves integrating with the right density of states for 



hyperboloids in D — 1 dimensions which is p2 

2 

^2(p) = TT-tanhTrp, Vs{p) ^ 



27r ' 27r2 

^^-^^^^ = 2n{D-?,) ^^-^^^^ = (A^^^v(D^^ ^\vUvW ' (3-24) 



2 


r(ip+f -i)|2 


(47r)^r(V) 


r(ip)|2 



We also need to define the regularized volumes of hyperbolic space in D — 1 dimensions. 
They are related to the volume of spheres 

^ ^ even 

^"^^ ^^"-^ Dodd ^^^^ 
2 

When D is even, we defined this regularized volume as minus the coefficient of log?]. 
When D is odd, we defined it to be the finite part after we extract the divergent terms. A 
derivation of these volume formulas is given in appendix B. Then the final expression for 
any dimension is 

Sintr = VHD-i,Teg dpV D-i{p) S{p, u) (3.26) 

^0 



with the expressions in ( p:25| ), {^^), (|3T9|) , (|3T7|) , ( [3:231) . We have defined ^reg as 



S =>S'intr log t] + ■ ■ ■ for D even 

(3.27) 

S =>S'intr + ■ ■ ■ for L> odd 

where the dots denote terms that are UV divergent or that go like powers of r] for small ry. 

3.6. Renyi Entropies 

We can also use the density matrix to compute the Renyi entropies, defined as: 

5, = -^logTVp^ g>0 (3.28) 
1 — q 



We first calculate the Renyi entropy associated to each SO{l, 3) quantum number. It 
is given by: 

S,{p, ^) = T^ log(l - l7pP) - log(l - |7pl"^) (3.29) 

1—0 1—0 
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Then, just like we did for the entanglement entropy (which corresponds to g — )■ 1), one 
integrates ( p.29|) with the density of states for D — 1 hyperboloids: 

/"OO 

-S'g.intr = Vfl-D-l^i-eg / dpV D-l{p) Sq{p, 1^) (3.30) 

With Sq^intr being the finite term in the entropy, for odd dimensions, and the term that 
multiplies log?7, for even dimensions. 

3. 7. Consistency checks: conformally coupled scalar and large mass limit 



As a consistency check of (|3.20 ), we analyze the cases of the conformally coupled 



scalar, and of masses much bigger than the Hubble scale. 
Conformally coupled scalar 

For the conformally coupled scalar in any dimensions we need to set the mass param- 
eter to u = 1/2. The entropy should be the same as that of fiat space. For a spherical 
entangling surface, the universal term is log euv I R for even dimensions, and is a finite 
number, g^ for odd dimensions |2y,^. The only difference here is that we are following 



a surface of constant comoving area, so its radius is given by i? = Rc/{Hr]). So, one sees 
that the term that goes like log 77, in even dimensions, has the exact same origin as the 
UV divergent one; in particular, we expect ce = ge for the four dimensional case, and go 
is the finite piece in the three dimensional case. 

Four dimensions: 

The entropy is given by ( p.22| ) 



This indeed coincides with the coefficient of the logarithm in the fiat space result pC 
Three dimensions: 
The entropy is given by: 



'S'intr =Vh2 e / f ^ tauh TT^g (p, ) = " / pdptanhnpS ( p, ^ 



'"^^ J (27r)2 ^ V 2 ; Jo V 2 

3C(3) log(2) 



(3.32) 



167r2 



This corresponds to half the value computed in |2^], because there a complex scalar is 
considered, and also matches to half the value of the Barnes functions in . 
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Conformally coupled scalar in other dimensions 

For even dimensions, 5'intr has been reported for dimensions up to (i = 14 in ||2^, and 
for odd dimensions, numerical values were reported up to d = 11 [^. Using ( p.26| ) we 
checked that the entropies agree for all the results in 12^,^ . 

Large mass limit 

Here we show the behavior of the entanglement entropy for very large mass, in three 
and four dimensions. The eigenvalues of the density matrix as a function of the S0(l,3) 



Casimir are given in terms of (|3.17| ). For large mass, there are basically two regimes, 
< p < \u\ and p > 

l7pl =|,-2.p (3.33) 
In this regime we can approximate I7I <^ 1 everywhere and the entropy per mode is 

^(p)~-|7pl'log|7pl' (3.34) 

Most of the contribution will come from the region p < up to 1/u corrections. 
This gives 



'-^^ ~ dpV{p)S{p) ~ (27rz.e-2-^) j^^ dpV{p) = | ^''^ ^ ^ ^ (3.35) 



H-D-i,reg 



which is accurate up to multiplicative factors of order (1 + C(l/z^)). 

3.8. Low mass range: l/2<z/<3/2 

In this low mass range the expansion of the field involves an extra mode besides the 
ones we discussed so far |[T^]. This is a mode with a special value of p. Namely p = i{u— 
This mode is necessary because all the other modes, which have real p, have wavefunctions 
whose leading asymptotics vanish on the S'^ equator of the future boundary. This mode 
has a different value for the Casimir (a different value of p) than all other modes, so it 
cannot be entangled with them. So we think that this mode does not contribute to the 
long range entanglement. It would be nice to verify this more explicitly. 

Note that we can analytically continue the answer we obtained for < 1/2 to larger 
values. We obtain an answer which has no obvious problems, so we suspect that this is the 
right answer for the entanglement entropy, even in this low mass range. The full result is 
plotted in fig. 2, and we find that for = 3/2, which is the massless scalar, we get exactly 
the same result as for a conformally coupled scalar. 
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4. Entanglement entropy from gravity duals. 

After studying free field theories in the previous section, we now consider strongly 
coupled field theories in de Sitter. We consider theories that have a gravity dual. Gauge 
gravity duality in de Sitter was studied in [^^,|2|,^||,|2|,^^,|33,^g|3^,^§7|,^, 
and references there in. When a field theory has a gravity dual, it was proposed in [rD| 
that the entanglement entropy is proportional to the area of a minimal surface that ends 
on the entangling surface at the AdS boundary. This formula has passed many consistency 



checks. It is certainly valid in simple cases such as spherical entangling surfaces [39|. Here 
we are considering a time dependent situation. It is then natural to use extremal surfaces 



but now in the full time dependent geometry . This extremality condition tells us how 
the surface moves in the time direction as it goes into the bulk. 

First, we study a CFT in de Sitter. This is a trivial case since de Sitter is conformally 
fiat, so we can go to a conformal frame that is not time dependent and obtain the answer 
ly,^. Nevertheless we will describe it in some detail because it is useful as a stepping 



stone for the non-conformal case. We then consider non-conformal field theories in some 
generality. We relegate to appendix C the discussion of a special case corresponding to 
a non-conformal field theory in four dimensions that comes from compactifying a five 
dimensional conformal field theory on a circle. 

4.1. Conformal field theories in de Sitter 




AdS, 



Fig. 3: The gravity dual of a CFT living on dS4. We slice AdS^ with dS'4 slices. 
Inside the horizon we have an FRW universe with i?^ slices. The minimal surface 
is an that lies on a constant global time surface. The red line represents the 
radial direction of this H^, and the S'^ shrinks smoothly at the tip. 
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As the field theory is defined in dS4, it is convenient to choose a dS^ shcing of AdS^. 
These slices cover only part of the spacetime, see fig. 3. They cover the region outside the 
hghtcone of a point in the bulk. The interior region of this lightcone can be viewed as an 
FRW cosmology with hyperbolic spatial slices. 

We then introduce the following coordinate systems: 
1. Embedding coordinates 



y^i - + Yl + ... + Yl = -1 

.2 _ jv2 jv2 , jv2 , , jv2 



(4.1) 



ds^ = -dYl^ - dY^ + rfy/ + ... + dYl 

2. dS/i^ and FRW coordinates 
2.1 dS slices 

YLi = cosh p, Yq = sinh p sinh r, Y^ = sinh p cosh rrii 

(4.2) 

ds'^ = dp^ + sinh p^—dr"^ + cosh T{da^ + cos^ adfl2)) 



2.2 FRW shces. We substitute p = ia and r = — + x in (|42 



y_i = coso", Yo = sin (J cosh X, Yi = sin a sinh xn^ 

(4.3) 

ds"^ = —da^ + sin^ o-{dx^ + sinh x{da^ + cos^ a(i02)) 
3. Global coordinates 

Y-i = cosh Po cos To, Yo = cosh p„ sin To, Yi = sinh pgnj 

(4.4) 

ds^ = dp^ — cosh^ PgdTg + sinh^ Pg{da^ + cos^ a(i02) 

As the entangling surface we choose the 5"^ at a = 0, at a large time tb and at p = cxd. 
In terms of global coordinates the surface lies at a constant r^, or at 

y 

sinhrs = tanr^s, Y4 = (4-5) 



Yli 
Its area is 

sinh^ pgdpg ~ 47r i-^ - ^ J (4.6) 

where pgc is the cutoff in the global coordinates. It is convenient to express this in terms 

of the radial coordinate in the dS slicing using sinhp^ = sinh p cosh r. In the large Pgc, Pc, 

Tb limit we find Pgc ^ Pc + '^b — log 2. Then (|4.6| ) becomes 

/g2pe+2Ts 1 \ / 1 1 \ 

A^4n ^-^^ -(p. + tb) j ~ 4. + 2 + ) j (4.7) 

We see that the coefficients of the two logarithmic terms are the same, as is expected in 
any CFT. Here euv = is the cutoff in the de Sitter frame and r] ~ e~^^ is de Sitter 
conformal time. 



17 



4-2. Non-conformal theories 

A simple way to get a non-conformal theory is to add a relevant perturbation to a 
conformal field theory. Let us first discuss the possible Euclidean geometries. Thus we 
consider theories on a sphere. In the interior we obtain a spherically symmetric metric and 
profile for the scalar field of the form 



(4. 



Some examples were discussed in pl| , p6| , ^ i. If the mass scale of the relevant perturbation 
is small compared to the inverse size of the sphere, the dual geometry will be a small 
deformation of Euclidean AdSo+i- Then we find that, at the origin, a = p + 0{p^), 
and the sphere shrinks smoothly. In this case we will say that we have the "ungapped" 
phase. For very large p we expect that log a oc p, if we have a CFT as the UV fixed point 
description. 



i 

FRW 

Crunch Hor 


dS / 

garn ^ 


Sac 


idle 


■9m 



Point 



Fig. 4: The typical shape for the scale factor for the gravity dual of a CFT per- 
turbed by a relevant operator in the "ungapped" phase. The region with negative 
corresponds to the FRW region. In that region, we see that = —a? reaches a 
maximum value, Om, and then contracts again into a big crunch. 



On the other hand, if the mass scale of the relevant perturbation is large compared to 
the inverse size of the sphere then the boundary sphere does not have to shrink when we 
go to the interior. For example, the space can end before we get to a = 0. This can happen 



^ We are interpreting the solutions of j42[j as explained in appendix A of ||12[j. This geometry 



also appears in decays of AdS space |42,42|. 
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in multiple ways. We could have an end of the world brane at a non-zero value of a. Or 
some extra dimension could shrink to zero at this position. This typically happens for the 
holographic duals of theories with a mass gap, especially if the mass gap is much bigger 
than H. We call this the "gapped" phase. See [l28| , ^^ , |33| , ^^ , |45[] for some examples. In 
principle, the same field theory could display both phases as we vary the mass parameter 
of the relevant perturbation. Then, there is a large phase transition between the two 
regime^. 

As we go to lorentzian signature, the ungapped case leads to a horizon, located at 
p = 0. The metric is smooth if a = p^-0{p^). The region behind this horizon is obtained by 
setting p = ia in ( [1.8| ) and (iO^ — )■ —ds\^ . This region looks like a Friedman-Robertson- 
Walker cosmology with hyperbolic spatial sections. 



ds^ = -da'^ + {a{a)fds]j^ , a{a) = -ia{ip) (4.9) 

If the scalar field is non-zero at p = we typically find that a singularity develops at a 
non-zero value of cr, with the scale factor growing from zero at a = and then decreasing 
again at the big crunch singularity. The scale factor then achieves a maximum somewhere 
in between, say at am- See fig. 4. 

We can choose global coordinates for dSn 

dsdSo = —dr"^ + cosh^ r(cos^ adQD-2 + da'^) (4-10) 

We pick the entangling surface to be the S^~^ at a = and some late time tb- We assume 
that the surface stays at a = as it goes into the bulk. In that case we simply need to 
find how r varies as a function of p as we go into the interior. We need to minimize the 
following action 

S = [ (acoshr)^-Vf^P^ - a'^dr^ (4.11) 

4GAr J 

The equations of motion simplify if we assume r is very large and we can approximate 
coshr ~ ^e"^. In that case the equations of motion give a first order equation for y = 

^ Since we are at finite volume we might not have a true phase transition. In de Sitter, thermal 
effects will mix the two phases. We will nevertheless restrict our attention to one of these phases 
at a time. 
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4-3. Non-conformal theories - gapped phase 

In the gapped phase, we can solve the equation for y. Inserting that back into the 
action will give an answer that will go like e*^^"^-*"^^ times some function which depends 
on the details of the solution. Thus, this produces just an area term. We can expand the 
action in powers of e"^"^ and obtain corrections to this answer. However, if the solution 
is such that the range of variation of r is finite in the the interior, then we do not expect 
that any of these corrections gives a logarithmic term (for even D) or a finite term (for 
odd D). Thus, in the gapped phase we get that 

^intr = (4.12) 



to leading order. The discussion is similar to the one in pq , ^ for a large entangling 
surface. 

4.4- Non conformal theories - ungapped phase 

In the ungapped phase, something more interesting occurs. The surface goes all the 
way to the horizon at p = 0. Up to that point the previous argument still applies and we 
expect no contributions to the interesting piece of the entropy from the region p > 0. 

When the surface goes into the FRW region note that the Sd-2 can shrink to zero at 
the origin of the hyperbolic slices. If we call p = ia and r = x "~ ^^/2, then we see that 
the metric of the full space has the form 

ds'^ = -da"^ + {a{a)f[dx^ + sinh^ ^(cos^ adnD-2 + da^)] (4.13) 

where d{a) = —ia{ia) is the analytic continuation of a{p). We expect that the surface 
extends up to x = where the S^~^ shrinks smoothly. Clearly this is what was happening 
in the conformal case discussed in the previous subsection. Thus, by continuity we expect 
that this also happens in this case. 

More explicitly, in this region we can write the action (|4.11|) , 



s.^/(..„.,)-J-(g)%a^ (4.14) 
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Fig. 5: The holographic setup for a non conformal field theory on de Sitter in the 
"ungapped" phase. We again have a region with dS4 slices and an FRW region 
with hyperbolic, H4, slices. The extremal surface that computes the entanglement 
entropy goes through the horizon into the FRW region. There it approaches the 
slice with maximum scale factor a = am- 

If we first set ^ = 0, we can extremize tlie area by sitting at cr^ where d = d^, which 
is the maximum value for d. We can then include small variations around this point. We 
find that we get exponentially increasing or decreasing solutions as we go away from am- 
Since the solution needs to join into a solution with a very large value of tb, we expect 
that it will start with a value of a at x = which is exponentially close to <Jm- Then the 
solution stays close to am up to x ~ tb and then it moves away and approaches cr ~ 0. 
Namely, we expect that for cr ~ the solution will behave as x = tb — log cr +rest, where 
the rest has an expansion in powers of e~^^^ . This then joins with the solution of the form 
r = — log p+rest in the p > region. The part of the solution which we denote as "rest" , 
has a simple expansion in powers of e~^'^^ , with the leading term being independent of tb- 
AU those terms are not expected to contribute to the interesting part of the entanglement 
entropy. The qualitative form of the solution can be found in figure fig. 5. 

The interesting part of the entanglement entropy comes from the region of the surface 
that sits near am- In this region the entropy behaves as 



Here we got an answer which basically goes like the volume of the hyperbolic slice H . 
This should be cutoff at some value x ~ tb. We have extracted the log term or the finite 
term, defined as the regularized volume. 




N 



intr 




reg 



(4.15) 
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Thus ( [4.15 ) gives the final expression for the entanglement entropy computed using 
the gravity dual. We see that the final expression is very simple. It depends only on the 
maximum value, am, of the scale factor in the FRW region. 

Using this holographic method, and finding the precise solution for the extremal sur- 
face one can also compute the coefficient C5 in (|2.6|) (or analogous terms in general dimen- 
sions). But we will not do that here. 

In appendix C we discuss a particular example in more detail. The results agree with 
the general discussion we had here. 



5. Discussion 

In this paper we have computed the entanglement entropy of some quantum field 
theories in de Sitter space. There are interesting features that are not present in the fiat 
space case. In fiat space, a massive theory does not lead to any long range entanglement. 
On the other hand, in de Sitter space particle creation gives rise to a long range contribution 
to the entanglement. This contribution is specific to de Sitter space and does not have a 
fiat space counterpart. We isolated this interesting part by considering a very large surface 
and focusing on the terms that were either logarithmic (for even dimensions) or constant 
(for odd dimensions) as we took the large limit. 

In the large area limit the computation can be done with relative ease thanks to a 
special S0(1,D-1) symmetry that arises as we take the entangling surface to the boundary 
of dSo- For a free field, this symmetry allowed us to separate the field modes so that 
the entanglement involves only two harmonic oscillator degrees of freedom at a time. So 
the density matrix factorizes into a product of density matrices for each pair of harmonic 
oscillators. The final expression for the entanglement entropy for a free field was given 
in ( |3.26| ). We checked that it reproduces the known answer for the case of a conformally 
coupled scalar. We also saw that in the large mass limit the entanglement goes as e~"^/^ 
which is due to the pair creation of massive particles. Since these pairs are rare, they do 
not produce much entanglement. 

We have also studied the entanglement entropy in theories that have gravity duals. 
The interesting contribution to the entropy only arises when the bulk dual has a horizon. 
Behind the horizon there is an FRW region with hyperbolic cross sections. The scale 
factor of these hyperbolic cross sections grows, has a maximum, and then decreases again. 
The entanglement entropy comes from a surface that sits within the hyperbolic slice at the 
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time of maximum expansion. This gives a simple formula for the holographic entanglement 
entropy ( [4.15|) . From the field theory point of view, it is an A^^ term. Thus, it comes from 
the long range entanglement of colored fields. It is particularly interesting that the long 
range entanglement comes from the FRW cosmological region behind the horizon. This 
suggests that this FRW cosmology is indeed somehow contained in the field theory on de 
Sitter space . More precisely, it is contained in colored modes that are correlated 

over superhorizon distances. 

In the gapped phase the order N"^ contribution to the long range entanglement entropy 
vanishes. We expect to have an order one contribution that comes from bulk excitations 
which can be viewed as color singlet massive excitations in the boundary theory. From 
such contributions we expect an order one answer which is qualitatively similar to what 
we found for free massive scalar fields above. 
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Appendix A. Bogoliubov coefficients 

Here we give the explicit form of the coefficients in ( p.lO|) . 



R ~ 



T{u + ip+ 1/2) ' ^ T{i^ + ip+ 1/2) 



T{u-ip + 1/2) ' "^"^ T{u-ip + 1/2) 



We also find 



r(,, J- 4- l)jp'^p+'i'^i' / 1 1 

^. ^ tP ^ 2'''^ I ie'^v+i-^'^ + 1 ie^p+in,._i 

4 sinh Trp \ j^ttp+TttIt^Ii j^^Trp+i-m^ _i / ju- 

r(';y — 7r) -I- / I I 

"''^ Asinhivp 



these were used to obtain ( |3.14| ). 
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(Al) 



(A.2) 



We define cr and cl via (|3.16|) and the state in ( p.l5|) . We demand that Ci^|\E') = 
7cJ^|\l/), cl|\E') = 7c^|\l/). Using (|3.16|) and denoting niRR = rriLL = P, = C, these two 



conditions become 



(up + V — 7f + [uC, — 'fvp — 'yu)b^j^ 








(A.3) 



(wC — 7W — 7fp)6^ + (wp + f — 7t'C)^l, 
which imply that each of the coefficients is zero. 

From the structure of ( |A.3| ), one sees that under the substitution u v v we 

have the same set of equations. If one tries to solve them together then ^ = hence this 
ratio must be real. One can show that this is indeed the case and 7 is given by ( p.l7|) . 



Appendix B. Regularized volume of the Hyperboloid 

Here we calculate the regularized volume of a hyperboloid in D — 1 dimensions. We 
have to consider the cases of D even and D odd separately. First, note that the volume is 
given by the integral: 



= VgD~2 

Now we expand the integrand: 

Vhd-1 1 



(ip(sinh p) 



D-2 



)D-2 



(^_l-^n^(D-2-2n)p 



(B.l) 



(B.2) 



But the integral of any exponential is given by: 



dp e 



ap 



1 

+ 

a 



0, a < 

divergent, a > 



(B.3) 



Now we treat even or odd dimensions separately. 

Even D: Here, the integrand of ( [B.2D contains a term independent of p in the sum- 
mation, which gives rise to the logarithm (a term linear in pc) . The term we are interested 
in corresponds to setting n = D/2 — 1: 



HD-i,reg 



(-1)^ fD-2 



D-2 \ D-2 



(-l)f7r^ {D-2y. (-l)f 
VsD-2 - /n-9x,o - 



n 



(B.4) 



2 / 2 

Odd D: Now, there is no constant term in the integrand of ( [B.2| ). Performing the 
summation in ( p.2| ), and using (|B.3|) , we get: 

D-2 



1 



2D-2 



E 

n=0 



D-2 



n 



i-iy 



D-2-2n 



D-2 



D-2 



-1 2 



(B.5) 
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A more direct way to relate the regularized volumes of hyperbolic space to volume of 
the corresponding spheres is by a shift of the integration contour. We change pc — > Pc + ^^r. 
This does not change the constant term, but we get an iix from the log term. We then 
shift the contour to run from p = along p = iO, with < ^ < tt and then from in to 
ZTT + Pc- The 6 integral gives the volume of a sphere and the new integral with Im{p) = n 
gives an answer which is either the same or minus the original integral. The fact that these 
regularized volumes are given by volume of spheres is related to the analytic continuation 
between AdS and dS wavefunctions p9|,[50[ . 



Appendix C. Entanglement entropy for conformal field theories on dS^ x 

Let us first discuss the gravity dual in the Euclidean case. The boundary is S"^ x S^. 
This boundary also arises when we consider a thermal configuration for the field theory 
on S"*. We will consider antiperiodic boundary conditions for the fermions along the S^. 
There are two solutions. One is AdS with time compactified on a circle. The other is the 
Schwarzschild AdS black hole. Depending on the size of the circle one or the other solution 
is favored [51,52|. Here we want to continue S"^ — > dS"^. An incomplete list of references 
where these geometries were explored is [p8i|30| , |3T| , ^|38| , |4^ , ^ . 

As a theory on dS'^ we have a scale set by the radius of the extra spatial circle. At 
large N we have a sharp phase transition. At finite we can have tunneling back and 
forth between these phases. Here we restrict attention to one of the phases, ignoring this 
tunneling. The Schwarzschild AdS solution looks basically like the gapped solutions we 
discussed in general above. Here, the 5"^, or the dS'^, never shrinks to zero. It can be viewed 
as a bubble of nothing. On the other hand, the periodically identified AdS^ solution gives 
the ungapped case, with the 5''* or dS^ shrinking, which leads to a hyperbolic FRW region 
behind the horizon. 

Gapped phase - Cigar geometry 
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Fig. 6: The gravity dual of a 5D CFT on dS4 x . The spacetime ends at r = r^, 
where the circle shrinks in a smooth fashion. We display an extremal surface going 
from tb to the interior. 



We now consider the cigar geometry. In the UV we expect to see the divergence 
structure to be that of a 5D CFT, but in the IR it should behave like a gapped 4D non- 
conformal theory. The metric is given by 

ds' = fd<p' + r'dsls. + ^, f = ^ + r'-^ (C.l) 
The period /3 of the 4> circle is given in terms of r^, the largest root of /(r/j) = 0, by 

Note that Pmax = tt/v^- This solution only exists for j3 < j3max- This geometry is shown 
in fig. 6. We consider an entangling surface which is an S"^ at a large value of tb- 
We need to consider the action 



A (X drr cosh ta/I — /r^r'^ (C.3) 



This problem was also discussed in [^. Since we are interested in large tb we can ap- 
proximate this by 



A 

-^approx 



j drr'^e^^^l - fr^ir')^ (C.4) 
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If the large tau approximation is valid throughout the solution then we see that the 
dependence on tb drops out from the equation and it only appears normalizing the 
action. In that case the full result is proportional to the area, e^^^, with no log- 
arithmic term. In the approximation ( |C.4| ), the equation of motion involves only r' 
and t". So we can define a new variable y = t' and the equation becomes first or- 
der. One can expand the equation for y and get that y has an expansion of the form 
y = [-2/(3r3) + 10/(27r^) - Am/{lAr^) + ■ ■ ■] a{l/r^ + ■ ■ ■) where a is an arbitrary 
coefficient representing the fact that we have one integration constant. 

This undetermined coefficient should be set by requiring that the solution is smooth 
at r = r/i. If one expands the equation around r = r^, assuming the solution has a power 
series expansion around rh, then we get that y should have a certain fixed value at rh and 
then all its powers are fixed around that point. Notice that if y = -t- y^(r — r^) -|- ■ ■ ■, 
implies that r is regular around that point, since (r — rh) oc where x is the proper 
distance from the tip. 

The full solution can be written as 



where y is independent of tb- 
A 



T = tb - I y{r)dr 



(C.5) 




rh 



Fig. 7: The regulated area A^eg is defined by Atotai = e^^^ Areg + Adiv 



At large r we get t — tb = + 0(l/r^) and the action (|C.4| ) evaluates to 



■^approx ^ 6 



2tb 



dr 



r^ + -+0 



.3 

c 

3" ■ 9 



~ e^^s ( + trc + finite 



1 — p^TB A 4- /I ,. 

(C.6) 
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We see that we get the kind of UV divergencies we expect in a five dimensional theory, as 
expected. 

These can be subtracted and we can compute the finite terms. These are plotted in 
fig. 7 as a function of rh- 

So far, we have computed the finite term that grows like the area. By expanding (|C.3| ) 



to the next order in the e~^^ expansion we can get the next term. The next term will give 
a constant value, independent of the area. In particular, it will not produce a logarithmic 
contribution. In other words, there will not be a contribution proportional to tb- 

In conclusion, in this phase, there is no logarithmic contribution to the entanglement 
entropy, at order 1/Gn or A^^. 

Ungapped phase - Crunching geometry 

Now the geometry is simply AdSg with an identification. This construction is de- 
scribed in detail in . The resulting geometry has a big crunch singularity where the 



radius of the spatial circle shrinks to zero. This geometry is sometimes called " topological 
black hole" , as a higher dimensional generalization of the BTZ solution in 3D gravity. 

It is more convenient to use a similar coordinate system as the one used to describe 
the cigar geometry in the previous case. The metric is given by ( |C.1|) , with m = 0. Those 



coordinates only cover the region outside of the lightcone at the origin, r = 0. To continue 
into the FRW region, one needs to use r = ia and r = —in/2 + % in (|C.1|) . 



The equation in the r > region is such that we can make the large r approximation 
and it thus reduces to a first order equation for y = r' = For small r, an analysis of 
the differential equation tells us that 

y~ 2r^H h6frH H j (C.7) 

where there is only one undetermined coefficient (or integration constant) which is b (it is 
really non- linear in b). This leads to a r which is 

r ~ -log r-r^H 'rc + b(r'^ + — —r^^ ) (C.8) 



where c is a new integration constant. We expect that the evolution from this near 
horizon region to infinity only gives a constant shift. In other words, we expect that 
c = TB+constant. This constant appears to depend on the value of b that is yet to be 
determined. We find that b should be positive in order to get a solution that goes to 
infinity and is non-singular. 
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We are now supposed to analytically continue into the FRW region. For that purpose 
we set r = ia and r = — z7r/2 + x- Thus the equation ( |C.8| ) goes into 



X ~ - log a - + ■ • • + c + 6 ( -a^ + H_[^^4 ^ . . . ) 



Then we are supposed to evolve the equation. It is convenient to change variables and 
write the Lagrangian in terms of <j{x) 

A~ j dx(T^siiv\v^x^-{(y'? + (^^0--(y^) (C.io) 

In this case, at x = we can set any starting point value for a(x = 0) and we have 
to impose that a' = 0. Then we get only one integration constant which is cr(0), as the 
second derivative is fixed by regularity of the solution to be (j"{0) = — cr(0)(3 — 4cr(0)^). 
We see that its sign depends on the starting value of ct(0). 

Since we want our critical surface to have a "large" constant value when we get to 
cr —7- as X —7- oo, we need to tune the value of a(0) so that it gives rise to this large 
constant. This can be obtained by tuning the coefficient in front of cr(0). This critical 
value of a(0) is easy to understand. It is a solution of the equations of motion with 
cr'{p) = (for a constant a), it is a saddle point for the solution, located at a = V3/2. 
If a is slightly bigger than the critical value, the minimal surface will collapse into the 
singularity, so we tune this value to be slightly less than the critical point. 

So, in conclusion, we see that the surface stays for a while at a ~ which is the 

critical point stated above. The value of the action ( |C1UD in this region is then 



3v^ p , . ,2 3^3 
/ ciY smh Y = 

16 io 16 



_Xc 
8 2 



(C.ll) 



Here Xc ^ t~b ^ \ogr] is the value of x at the transition region. Thus, we find that the 
interesting contribution to the entanglement entropy is coming from the FRW region. 

The transition region and the solution all the way to the AdS boundary is expected 
to be universal and its action is not expected to contribute further logarithmic terms. 

In conclusion, the logarithmic term gives 
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Here we repeat this computation in another coordinate system which is non-singular 



at the horizon. We use Kruskal-hke coordinates P^ , ^|30| . It also makes the numerical 
analysis much simpler. In terms of embedding coordinates for AdS^: 

- y2 _ y2 ^ y2 ^ ^ y2 _ 

(C.13) 

ds^ = -dY^^ - dY^ + dY^ + ... + dYi 
The Kruskal coordinates are given by: 

Y-i = - ^cosh(/), Y5 = - jsmh(/), yo,...,4 = 5-, 2/ = -2/o + 2/i + ••• + 2/4 

1 — 1 — 1 — y 

(C.14) 

In these coordinates, the ^5" region corresponds to < < 1 and the FRW region 
to —l<y^< 0, with the singularity located at = —1. The AdS boundary is at 
y'^ = 1. We can relate the pair (r, r) and (x,c7"), connected by the analytic continuation 
(r = zcr, r = — z7r/2 + x) to (y^,yo) by the formulas: 



r = , smh r 



^ (C.15) 
CT = — ^, coshx- 



The functional gets simplified to: 




2^2 



^= / \/ ^li^yl)^ + ?/^) [(^(j/^))^ + ^yodyodiy^) - 4y^{dyoy] (C.16) 



If one looks for the saddle point described in the FRW coordinates, then one obtains 
y^ = —1/3. The same situation can be described in simple fashion in terms of these 
coordinates. So ~ —1/3 for a large range of yo, and it crosses to the (i5'-sliced region 
at a time of order r^, in dS coordinates, so part of the surface just gives the volume of an 



if 3 



as m 



9 

16 



cosh Xc 
V3 



3yo - 1 dyo 



3^3 



Xc 

T 



(C.17) 



Some plots for the minimal surfaces are shown in fig. 8. 
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Fig. 8: We plot here the value of versus j/o- For small yo the solution starts 
very close to the surface of maximum expansion at = —1/3, stays there for a 
while and then they go into the AdS boundary at y'^ = 1. The closer yo is to 
the saddle point y^ = —1/3, the longer the solution will stay on this slice, giving 
a contribution that goes like the volume of an . Then, at a time yo of the 
order of the time the surface reaches the boundary, it exits the FRW region. The 
interesting (logarithmic) contribution to the entropy is coming from the volume of 
the surface along the FRW slice at = —1/3. 
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